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Abstract. We present limits on transit timing variations and secondary eclipse depth variations 
at 8 microns with the Spitzer Space Telescope IRAC camera. Due to the weak limb darkening 

Cu . in the infrared and uninterrupted observing, Spitzer provides the highest accuracy transit times 
for this bright system, in principle providing sensitivity to secondary planets of Mars mass 
in resonant orbits. Finally, the transit data provides tighter constraints on the wavelength- 

^. ■ dependent atmospheric absorption by the planet. 

JjL ' Keywords, stars: planetary systems 

o 
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The extremely precise 33 hour phase function measurement of HD 189733b Knutson 

et al. (2007) observed with the 8 micron IRAC camera on the Spitzer Space Telescope 

yielded the most precise times of transit and secondary eclipse for any extrasolar planet, 

6 and 24 seconds, respectively. This led us to propose to observe an additional six transits 

i* . and eclipses of this system over time with the goal of measuring: 

• precise transit-timing (Agol et al. 2005, Holman & Murray 2005) to search for the 
presence of resonant (or near-resonant) terrestrial-mass planets captured by migration 
(e.g. Mandell, Raymond & Sigurdsson 2007) or on longer timescales precession of an 
eccentric orbit (Miralda-Escude 2002, Heyl & Gladman 2007; also Fabrycky and Wolfe, 
these proceedings); 

• variations in the depth of the secondary eclipses with time which might be caused 
by large-scale variable atmospheric features (Rauscher et al. 2007; also Showmand and 
Dobbs-Dixon, these proceedings); 

• precise transit depth for comparison with atmospheric-absorption models to contrain 
the molecular composition (e.g. Tinetti et al. 2007; also Tinetti et al., Fortney et al., and 
Hubeny et al., these proceedings); 

• improved system parameters for better characterization of the planet, host-star, and 
orbit properties (Winn et al. (2007); also Winn, these proceedings). 

For librating planets in a low-order mean motion resonance, the times of transit vary 
with an amplitude: 

r . D / Mpert \ _„ . ( Ptrans\ ( M pert \ ( Mj \ . . 

^ * Ptr - \M^nS) * 3 mm UdiyJ VW) U^J ' (L1) 
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and libration period of 

where M ffi , Mj, M* are masses of the Earth, Jupiter, and host star, Ptrans is the period of 
the transiting planet, and the numbers have been estimated for the libration amplitude 
for planets starting on circular orbits with exact commensurability (Agol et al. 2005); 
the actual value depends on the libration amplitude. This timescale requires observa- 
tions separated by months with k, seconds precision in timing, and in principle could be 
sensitive to sensitive to Mars mass planets. 

2. Data reduction summary 

So far Spitzer has observed four transits and four secondary eclipses of HD 189733 for 
this program with 44,000 exposures of 0.4 second each over 5 hours each. Additionally 
we re-analyzed the data from Knutson et al. (2007). We utilized the IRAC 8 microns 
as it has been demonstrated to be the most stable IRAC band. Due to the brightness 
of the host star we made the observations in sub-array mode. We carried out aperture 
photometry with a 3.5 pixel radius. Due to the small limb-darkening, stable instrument 
(thanks to the Earth-trailing orbit of Spitzer which leads to stable thermal properties and 
no occultation of targets by the Earth, as occurs with HST) we obtained 0.5% precision 
per 0.4 second exposure. In the shot-noise limit we expect a 3-second precision for transit 
times. 

Figure [1] shows a gallery of the transits and secondary eclipses obtained. The strong 
"ramp" which causes the flux to change by about 1% is a well known feature of the IRAC 
camera when observing bright sources. We fit the ramp with the function clq + a\(t — 
ts) + a 2 hi (t — t s ), where t s is the start of the observation and do, 1,2 are constants. This 
appears to remove any trace of the ramp in the corrected data. We assign each data point 
an error bar which equals the scatter in residuals of the data outside of transit/eclipse. 

3. Transit/eclipse depths and limb darkening 

In two initial fits we include one ephemeris for the transits, one ephemeris for the 
secondary eclipses, and required all other parameters describing the transits and eclipses 
to be the same for all 10 transits/eclipses to determine their mean values. We performed 
one of these fits with stellar limb-darkening set to zero, while for a second fit we allowed 
the linear coefficient of the stellar limb-darkening to freely vary (in all cases we assume 
the planet to be uniform in surface brightness). We allowed the ramp parameters to vary 
independently for each transit/eclipse. A sum of the five transits corrected for the ramp 
is shown in Figure EHa) without limb darkening, and Figure [3th ) with limb darkening. 
The five secondary eclipses are shown in Figure [3J 

Without limb darkening we find a best-fit x 2 = 232186.8 for 231508 degrees of free- 
dom (231516 data points with 8 model parameters). With linear limb darkening the % 2 
improves by A\ 2 = 140 for 231507 degrees of freedom, with a best-fit limb darkening 
parameter of U\ — 0.110 ± 0.015; thus limb-darkening is detected at 12cr - this can be 
seen by eye by comparing Figures [3J a) and[3Jb). This best-fit value for limb-darkening 
gives a limb-darkening profile which is very close to that predicted by a Kurucz model 
with parameters appropriate for HD 189733a (effective temperature T e // = 5000 K 
and surface gravity log[g(cm/s 2 )] = 4.5). We find a best-fit planet-star radius ratio of 
Rp/R* = 0.1558 ± 0.0002, which translates into a best-fit area ratio of 2.427 ± 0.005%, 
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while the best-fit eclipse depth is 0.347 ± 0.005% - note that this is a 72ct detection of a 
secondary eclipse! The errors on these parameters (and throughout the rest of the paper) 
are computed by generating one hundred synthetic data sets from the best fit light curve 
(for each model) added to the residuals shifted by a random number of data points, 
re-fitting the model to each synthetic data set, and then taking the standard deviation 
of the 100 synthetic parameter sets to determine each parameter's error. 

4. Ephemeris and Transit Timing Variations 

We have measured the best-fit ephemeris to our data, separately for the transits and 
eclipses. We find a transit ephemeris of T = 2454279.436741 ± 0.000023 HJD and P = 
2.21857503 ± 0.00000037 days, while for the eclipses we find T = 2454279.437407 ± 
0.000130 HJD and P = 2.21857306 ± 0.00000209 days, assuming that they occur exactly 
one half period after primary transit. The periods are consistent within 1 — a, while the 
central eclipse times differ by 57±11 seconds, indicating that the secondary eclipses occur 
later. About 31 seconds of this difference can be accounted for by the light travel time 
across the system, while the other 26 seconds is likely due to a slight orbital eccentricity. 

We have performed a second set of fits allowing the central times of transit/secondary 
eclipse to vary. The transit timing variations are plotted in Figure [4] for the transits 
and secondary eclipses, compared to the best-fit transit ephemeris (and assuming the 
secondary eclipse is offset by 1/2 an orbit). We find no significant deviations from a 
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Figure 1. Atlas of transits and secondary eclipses obtained at 8 microns with Spitzer. 
Horizontal axis is time in units of days; vertical axes are photon counts in arbitrary units. 
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Figure 2. (a) Binned transits with best-fit model (solid red line) assuming no limb darkening 
for star. Data include error bars, (b) Transits fit with best-fit model with limb darkening (red 
solid line). 

uniform period by more than 5 seconds for either the transits or secondary eclipses; once 
the final four data sets are obtained we will carry out a more detailed analysis. Poster 
79 by Miller-Ricci and poster 87 by Neale Gibson at this IAU symposium also present 
new transit-timing data for this system. 



5. Eclipse Depth Variation 

In a third set of fits we held the ephemerides fixed, but allowed the depths of secondary 
eclipse to vary. Figure [5] shows the variations in the five observed eclipse depths in units 
of the flux of the star. We find that the fractional variations in the eclipse depth are 
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Figure 3. Binned data from 5 secondary eciipses with best-fit model (red solid line). 
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smaller than about 10%; this is about the level of variation predicted by Rauscher et 
al. (2007). A fit to the eclipse depths assuming a constant depth gives a \ 2 — 6 for 4 
degrees of freedom (5 eclipses minus one free-parameter, the mean eclipse depth). Once 
the final two secondary eclipse observations are obtained (in June and July 2008), a 
more comprehensive analysis will be carried out of the limits on secondary eclipse depth 
variation. 



6. Transit spectroscopy 

The depth of primary transit we have obtained contains five times as much data as 
that in Knutson et al. (2007), providing more precise constraints on the spectral energy 
distribution in the mid-infrared. We have combined this data point with transit depths 
measured by Pont et al. (2008) in the optical with HST, Swain et al. (2008) in the 
near-infrared with HST, Beaulieu et al. (2008) at 3.6 and 5.8 microns with IRAC, and 
Knutson et al. (2008) at 24 microns with MIPS (also Knutson et al., these proceedings). 
The data are plotted in Figure [6] along with the best-fitting model of Tinetti et al. (2007) 
which does not have enough methane to fit the near-infrared data. We find that our 8 
micron data point lies above the Tinetti model; as the 8 micron band covers a region 
in which methane absorption is quite strong (~ 10 20 cm 2 per molecule), the model may 
be brought back into agreement with the data with a higher abundance of methane, as 
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Figure 4. Observed minus calculated (O-C) transit times for the five primary transits (black 
data points with small error bars) and the five secondary eclipses (black data points with larger 
error bars), compared to the ephemeris measured by fitting just the primary transits (red solid 
line) with uncertainties (red dotted lines). The ephemeris from fitting the secondary eclipses 
is shown as a thick blue solid line, with dashed blue lines showing the uncertainty in slope. 
Diamond green data points are transits from Pont et al. (2007) measured with Hubble Space 
Telescope with similarly sized error bars to our data. 
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is already required by the near-infrared data (Swain et al. 2008). It also clear that the 
model is a poor fit to the optical data; the fit may be improved by including Rayleigh 
scattering (Lecavelier des Etanges et al. 2008, and these proceedings). 



7. Conclusions 

The five transits we have observed are consistent with no transit timing variations 
greater than 5 seconds. The five secondary eclipses are consistent with no eclipse depth 
variations at greater than the 10% level. From just an analysis of the lightcurves we can 
place constraints on the eccentricity and longitude of periastron: |ecoscj| = 0.0002 ± 
0.0001 and |e sinw| = 0.015 ± 0.012; similar to the results of Winn et al. (2007). We are 
collecting two more transits and two more secondary eclipses in Summer 2008; thus the 
results presented here are preliminary and require further more careful analysis. 
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Figure 5. Changes in eclipse depths; cyan line is mean of 5 eclipses (0.347%). 
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Figure 6. Compilation of data on transit depths of HD 189733. Solid (red) curve is a binned 
version of a model from Tinetti et al. (2007) ; horizontal lines (black) are the mean of this model 
over the bandwidths of each of the crosses (blue), which indicate 1-a error bars. 



